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Abstract 


Maneuvering characteristics of surface combatants in the United States Navy are often 
ignored durmg the design process. Key maneuvering parameters such as tactical diameter 
and turning rate are determined during sea trials after the ship enters service. Iu the “Navy 
After Next”. the study of maneuvering of surface combatants will become increasingly more 
important in efforts to reduce the number of personnel required to operate the ship and 
thus reduce hfe cycle costs. This thesis attempts to address this issue. 

The thesis presents an Extended Kahnan Filtering (EKF) algorithm to estimate the 
linear damping hydrodynamic coefficients for an Arleigh Burke Class Destroyer. Actual 
data is generated by conducting maneuvers (witl: a nonlinear model of the ship developed 
In a separate study) where nonlinear effects are small. The EKF then uses that data 
to estunate the lull stability coefficients (Y,. Ny. Y;, and N,) on-line in real time. The 
cocfficieut values deternmned by the EKF are then used in a simulation model and. the 
results are compared to the actual trajectories. Despite the nonlinearities present in the 
actual data, the EKF provides coefficient vaiues that reproduce trajectories with only 15% 
OIOI: 

The linear coefficients are then used to develop simple controllers to automate maneu- 
vering for the actual ship. The parameters determined by the EKF are used to derive a 
linear tine invariant (LTI) model of the ship. This LTI model then serves as the basis 
for model-based compensator designs to automatically control ship maneuvers. The first 
controller 1s an autopilot to regulate the ship’s heading and the second is a regulator that 
ensures the ship remams on its intended track. The performance of the compensators 1s 
then evaluated by simulating the performance of the LTI controllers on the nonlmear plant. 


Thesis Supervisor: Michael S. Triantafyllou 
Title: Professor of Hydrodynamics and Ocean Engineering 
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Chapter 1 


Introduction 


1.1 Background 


In the present naval ship design environment, reduced manning represents the cornerstone 
of any design. Mannimg reductions significantly reduce the life cycle cost of the project. 
which is the main goal given today’s shrinking defense budgets. Reduced manning also 
lunplies fewer sailors being placed in harm’s way during battle as well as enhanced quality 
of life. Reducing manning. however. affects nearly all aspects of a ship: maintenance. 
hrefighting aud damage control. maneuvering. etc. The solution to this problem requires 
a major paradigm shift from the Navy s current doctrine, culture. tradition. and traming 
procedures. 

Naval ships must operate in a multitude of different environments and perfornı well m 
cach. Maneuvering of the ship in each of these environments currently requires numerous 
watchstanders on the bridge and in the Combat Information Center (CIC) during normal 
operations. Of these numerous watchstanders. one must be the helmsman. the individual 
respousible for manipulating the rudder and engines to keep the ship on its desired course 
and speed. With the advances in modern control theory and computing power, the functions 
performed by many of these watchstanders may potentially be automated. This represents 
a major change in current. Navy operational doctrine, however. 

Automation of these functions requires some method of automatic control. The evolu- 
tion of control can be broken into two distinct periods. The period prior to 1957 can be 
cousidered the “Classical” period and the period from 1957 to the present can be considered 


the “Modern” period [33]. Classical control theory deals mainly with single-imput /single- 





output (SISO) linear time-invariant (LTI) systems. whereas modern control theory expands 
that capability to deal with multiple-input/multiple-output (MIMO) systems. 

Automated ship maneuvering has been attempted since the inveution of the gyroscope 
by E. A. Sperry in 1910. Sperry applied his gyroscope to the stabilization and steering 
of ships and later aircraft [33] using classical control methodology. Maneuvering of naval 
surface combatants. however. is almost completely neglected in the design process. Naval 
architects are content to accept the ship’s maneuvering capabilities as determined by full- 
scale sea trials after the ship has been built. The manning goal of the USN's 21% century 
surface combatant (DD-21) is 95 personnel. Automation of maneuvering must surely be 


emploved in the desien if this goal is to be realized 


1.2 Objectives 


The main objective of this work is to apply a system identification technique known as 
the State Augmented Extended Kalman Filter (SAEKF) to identify the hydrodynamic 
coefficients of a model of the DDG-51 on-line in real time from noisy measurement. data. 
The data obtamed from the identification may then be used to develop simple. controllers 
that can be tested during the design process for implementation on the full-scale ship. 
It is worth noting that the simulations performed in this study may not match full-scale 
DDG-51 trial data. The model used in this study to generate the ship maneuvering data 
was developed in [43]. No attempt has been made to validate this model against. full-scale 
manenvering data. The goal is to properly identify coefficients in an assumed form of the 
maneuvering equations and attempt to reproduce the states and trajectories produced by 
the model from [43]. In terms of developing control laws. the most important factor is the 
values of the systein states. rather than exact values of the hydrodyuamic coefficients for a 
particular ship. In fact, linear control systems often perform quite well despite errors as large 
as 40% in the states. Therefore, if the identified coefficients produce values of the system 
states that are close to those generated by the model in [43], that form of the simulation 
equations can be assumed to be accurate and may then be used to develop control laws to 
automate the slip's maneuvering. 

In order to achieve the goal, a form of the simulation equations for ship maneuvering 


must be assumed. This form is then used to develop an SAEKF algonthin that accepts noisy 
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measurement data from ship maneuvers generated by the model in [43] and estimates the 
value of the hydrodynamic coefficients in an attempt to make the output of the simulation 
model match the actual measured data. The form of the model is based upon physical 
principles aud contains the salient terms required to describe the coupled surge. sway, and 
vaw motions of tle ship. Dimensionless quantities are employed throughout the process to 
maintain generality and ensure numerical stabihty of the algorithm. 

The identification portion of the study limits its focus to forward motions in deep, cali 
seas with no current. The controller stmulations. however, relax this condition and use 


curreuts aud waves as disturbances to test the performance of the controller designs. 


13 Contributions 


This work describes and evaluates a process by which coefficients in maneuvering equa- 
tions of motion can be identified on-line in real time. From these identified parameters. 
control laws for automating the ship's maneuvering in real time m the ship’s operational 
enviromuent can be developed. This method can be utilized early in the design process 
on scale models to deternime the effectiveness of the automatic control systems aud deter- 
mine changes that need to be made. This work can be extended to apply to future naval 
combatant ship designs such as DD-21. 

This work introduces a method by which naval architects can address maneuverimg 
characteristics early in the desigu process aud develop automatic controllers to reduce the 
number of human interfaces required to maneuver the ship. The method is written en- 
tively in the MATLAB computing environment. This work will also help to provide Ocean 
Engineering graduate students at MIT with practical experience in guidance and control 
of ocean vehicles through its potential future implementation on a scale model. Further, 
the simulations developed in the thesis may be used as the basis for future student design 


projects i slip maneuverimg and control. 


1.4 Outline 


Chapter two presents a truncated derivation of the governing equations for the SAEKF 
algorithm. It then presents au application of the identification process to a very large 


crude carrier (VLCC) Esso Osaka. This is an application of the significant work couducted 
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bv Professor Martin Abkowitz of MIT in the late 1970's and early 1980's. This model 
is developed aud used for verification of filter operation throughout the reniainder of the 
thesis. 

Chapter three presents a brief description of the United States Navy's Arleigh Burke 
Class (DDG-51) Destroyer. It then develops the governing equations of ship motion in tlie 
horizontal plane. The chapter concludes by developing the assumed form of the equations 
of motion to be employed in the SAEKF. 

Chapter four develops tlie linear equations of ship maneuvering motion. It describes 
the results of the initial identification efforts and problems encountered. The chapter fur- 
thor describes bias and divergence in the Kalman Filter and some of the causes of tliese 
phenomena. The methods of addressing the problems encountered are then described. The 
chapter concludes with the results of the successful identification of the linear daniping 
hydrodynamic coefficients for tlie DDG-51. 

Chapter five briefly describes the theory of controller design by loopshaping. Two sin1ple 
coutrollers are then designed using loopshaping and the linear models produced through the 
system identification efforts. The first controller is an autopilot designed to maintain the 
ships actual heading about a slowly-varying reference value. Simulations of the controller 
are performed on the nonlinear plant to determine its performance. Cross-track errors 
evidenced i the autopilot design then lead to the design of a siniple track-keeping controller. 
Simulations are performed oun the noulinear plaut for tlis desigu as well. 

Chapter six summarizes the work accomplished in the thesis and gives recommendations 


for future studies. 
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Chapter 2 


System Identification and the State 
Augmented Extended Kalman 
Filter 


This chapter introduces the use of the Extended Kalman Filter (EKF) in system identifica- 
tion. The govermng equations of the EKF appear along with a brief description of the noise 
processes. The chapter concludes with an example application of the EKF to the identifi- 
cation of unknown hydrodynamic coefficients for the very large crude carrier (VLCC) Esso 


Osaka. 


2.1 System Identification 


System Identification (SI) is the process of developing or improving a mathematical rep- 
resentation of a physical system using experimental data. This process usually assumes a 
form of a model for a physical system and adjusts the unknown paranıeters ın that ınodel 
to fit physical data collected from the system in question. The State Augmented Extended 
Kalınan Filter (SAEKF) illustrates one technique for performing SI. The system state vec- 
tor is augmented with the system’s unknown parameters and estimated as data is collected 
from the physical system. 

Many methods exist for performing SI on mathematical models of physical systems. 


Several of these methods have been successfully applied to the ship maneuvering problem. 
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Maximum Likelihood Parameter Estimation was successfully applied in [41] to identify 
linear coefficients. Reference [13] covers several other methods such as Indirect Model 
Reference Adaptive System. Continuous Least Squares Estimation, and Recursive Least 
Squares Estimation. Abkowitz [3] successfully applied the SAEKF inethod to tanker ships 
i the early 1980°s. The SAEKF method provides a means to update a system model on-line 
i real tune. The success of the work presented in [3], coupled with the real time estimation 


capability. forms the basis for the choice of the SAEKF method in this study. 


2.2. The Extended Kalman Filter 


In the most general case. nonhnear systems and subsequent sets of discrete measurements 
of those systems can be described by a set of nonlinear. stochastic differential equations of 


the following fornı !: 


where: 

w(t) represents the svstein state vector. 

f(.c(t).t) represents the system description matrix. 

w(t) represents a zero-mean Gaussian sequence with covariance matrix Q(t), 
zy, represents the measurement vector. 

hy (.r(t,)) represents tlle measurement description matrix, and 


ey, represents a zero-mean Gaussian sequence with covariance matrix Rx. 


fu the ship maneuvering problem, system measurements typically cousist of direct mea- 
surements of the state variables. Thus. the matrix hy,(z(t,)) reduces to a constaut identity 
matrix. For example. ships often utilize a Global Positioning System (GPS) to directly 
measure the ship's position and speed. The ship also has a gyrocompass to directly mea- 
sure the instantaneous heading angle. The direct state measurements require no additional 
calculations to determine the value of the state. Therefore, the system measurement matrix 


'The derivations that follow appear in more detail in [19]. 
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reduces to a constant function equal to the identity matrix.” Keeping this fact in mind. 


equation 2.2 reduces to the following: 


Zk = hkz(tk) + vy 23] 


The equation of motion 2.1 aud the measurement equation 2.3 now govern the dynamics 
of the entire system. The EKF method seeks the minimum variance estimate of z(t) as a 
funetion of time and the measurement data accumulated up to time t. The sequel presents 
an abbreviated derivation of the theory behind the EKF algorithm. Reference [19] presents 


a more detailed derivation of the interinediate steps. 


2.2.1 Process and Sensor Noise 


The noise sequences w aud v, in equations 2.1 and 2.3 represent the process aud sensor 
noise. respectively. Their properties aud associated covariance matrices have yet to be 
addressed. 

In the ship maneuvering problem. the process noise represents the uncertainty in the 
assumed form of the model as well as uncertainty in predicting external disturbances. For 
example. the sea mav not always exhibit a wave spectrum exactly consistent with that 
predicted by statistical data. In calm seas. however, this uncertainty is removed and only 
the uncertainty in the niodel forin remains. Thus, the process noise covariance may be 
assumed to be constant. This is not a limitation. however, because even in a high sea state. 
the additional uncertaity due to the wave excitation will vary slowly over tune and can 
be considered piecewise constant [32]. Furthermore. because the process noise represents 
uncertainty. it is assumed that the state vector and process noise are independent and 
uncorrelated random variables. 

The sensor noise represents the uncertainty in the measurements. For example, GPS 
systenis available in the commercial market provide position information to within +/- 100 
feet in some instances, but not the exact location. The measurement model must address 
this uncertainty. Therefore. in the sequel. the following assumptions hold for the process 


and sensor noise: 


> = . . . . . . 
“If all state variables are not measured, the matrix will consist of ones and zeros of sufficient size to 
perforin the required linear algebraic operations. 
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lL. White noise processes.” 


e The ship motions are slow compared to the dynamics of waves. structural vibrations. 
etc. Thus. the tinie constant of the process noise is much, much faster than the time 
constant of the ship motions. 

e The observation interval must be long compared to the correlation time of the sensor 
noise. Typical shipboard sensors exhibit the capability of sampling rates as high as 100 
hertz. Therefore. a sampling interval of not less than oue second meets this condition 


and is more than adequate to fully capture the slow ship motion dynamics. 


2. Independent. uncorrelated. zero-mean, Gaussian random variables (denoted N(0.Q)). 
7 I 


Thus. E futte] E ul)” | = (0 for all k and t. 
3. Covariances are constant or piecewise constant varying slowly witl time. 
These assumptions are quite important in the derivations that follow. Reference [19] pro- 
vides a more detailed discussion of the noise processes. 
2.2.2 State Estimate and Error Covariance Propagation 


The EKF method seeks to ininimize the error of the estimate in some statistical sense as a 
function of time. Thus. the error obviously depends on time. Define the error, :c(t), and its 


associated error covariance matrix. P(t) as follows: 


z(t) = Z(t) — alt) (2.4) 
P(t) = E[g(t)z(t)"] (2.5) 


where T(t) denotes the minimum variance estimate determined by the EKF. and Ef e ] 
denotes the expectation operator. From the definitions in equations 2.1, 2.4. and 2.5. it can 
be shown [19] that equations 2.6 and 2.7 govern the state estimate and its associated error 
covariance propagation. Equation 2.7 omits the time dependence of f and z for notational 


“Nature does not exhibit white noise processes. Thus, the white noise assumption is only valid subject 
to the conditions presented. 
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Convelence. 


(26) 


[> * 


(t) = F(x(t),t) 
P(t) = xf? - ET Y faT — Fat + QUe) (2.7) 


| 


IS) 


lu order to solve these equations, f must be known at each instant. The value of f, however. 


depends upon the entire probability density function of x at each instant of time as shown 


by the following relation: 


ee / i / A E, (2.8) 


It 


where p(x.f) denotes the probability density function of x. Practical algoritlms require 
methods of computing z and P without knowing p(x,t). One method of accomplishing this 
is by expanding f in a Taylor series about the current estimate, x. as follows (assuming the 
required derivatives exist): 


(2.9) 


Truncating equation 2.9 to include only the first order terms. taking the expected value. 
and substituting ito equations 2.6 and 2.7. yields the following approximate expressions 


for the propagation of the estimated state vector and its associated error covariance: 


Belt) (2.10) 
P(t) = F(&.t)P(t) + P(t)F’ (&.t) + Q(t) (2.11) 
where F(£.t) is the matrix whose ijt? element is defined as follows: 
Oa 
filz.t) a 


i 17 i 3 x; noe 


2.2.3 State Estimate Update, Error Covariance Update, and the Kalman 


Gain Matrix 


The filtering equations 2.10 and 2.11 propagate the state estimate and error covariance 


over the time period between measurements (t,_; < t < tk). Therefore, to complete the 
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filter aud update the state estimate for the next time step, the actual state measurement 
must be taken into account. Because the EKF algorithni produces a mininium variance 
estimate of the state vector, there will likely be a difference between the propagated state 
estimate (denoted £,(—) +) and the actual measurement (denoted z,). Thus. the updated 
state estimate (denoted 7,(+)) can be computed by a linear combination of the propagated 
state estiinate ald the difference between the actual measurement and propagated state 
estimate. The scaling factor in this linear combination is known as the Kalman Gain 


matrix. Ky. Equation 2.13 illustrates this relation. 
£y(+) = £,{(—) + Ky [z, - helz,(-))] (2518) 


Measurements affect the error covariance in a manner similar to the state estimate 
update. Therefore. denote the propagated error covariance as P,(—) and the updated 
error eovariauce as Py, (+). The optimum Kalman Gain matrix, Kk. minimizes the error 
covariance update. P,(+). Expressing P;,(+) as a function of of K, yields the following 


expression to be minimized: 
Pi (+) = Px(—) + Ky E| [zy — hy (a,)]2,(-)* (2.14) 


Like the state estimate update, T¿(+). equation 2.14 is a linear combination of the prop- 
agated error covariance and tlie estimation error. Mininlizing equation 2.14 yields the 


following expression for the optimum Kalman Gain matrix (see [19] for details): 


r 


2i 
MS AS 2, — fi ea E las == hi (ea) | + Rif a) 


A 


Note that hy, depends upon the entire probability density function of x(t) similar to f in 
2.2.2 Einploving the method in 2.2.2 and expanding h, in a Taylor series about. the 


current propagated state estimate. z,(—). yields the following expression: 


A A Ll) es (2.16) 


'r,(—) and P,(—) represent the solutions to equations 2.10 and 2.11 on the interval tp- < t < tp. 
“See equation 2.8. 
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where: 





d Phy, (x) 
Haad a Haa (2.17) 
Ox r=2,(—) 
Thus. similar to F. Hy is a matrix whose ij element is given by the following expression: 


(Hi); = (Zs) 





Ly =Z,(—) 


Truncating equation 2.16 after two terms, taking the expected value and substituting into 


equations 2.13. 2.14, and 2.15 yields the Extended Kalman Filter update equations: 


ee ee h(a —))| (219) 
= 

Ki =P JH (BA) E NP Ha (&(—)) + Re (2.20) 

BES E EEE PR) (22210) 


Figures 2-1 and 2-2 illustrate the process with a flowchart and timing diagrain. Table 2.1 
presents a summary of the Extended Kalman Filter algorithm. Each block in the flowchart 
represents a section from Table 2.1. Figures 2-1 and 2-2 should be used in conjunction with 
Table 2.1 to illustrate the flow of the algorithm and the governing equations for each stage 


of the process. 
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Figure 2-1: EKF Algorithm Flowchart 
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Figure 2-2: EKF Timing Diagram Adapted from Figure 4.2-1 of [19] 
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Table 2.1: Summary of the Extended Kalman Filter Equations 
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State Estimate 7 = f(z.t) (Plant Simulation Model) 


Propagation 





Error Covariance P(t) = F(%. t)P(t) + P(t)F? (z. t) + Q(t) 
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Adapted from Table 6.1-1 of [19] 
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2.3 System Identification of the Esso Osaka 


No a priori information on the parameters of the DDG-51 was available in this study. Thus. 
if the identification process produced erroneous values of the coefficients, there would be 
uo definitive method of isolating the problem to the model or the filter. In an effort to 
clisure proper operation of the filter prior to use in identifying the DDG-51 hydrodynamic 
coefficieuts. a linear model of the Esso Osaka was developed using coefficient values provided 
al. 

The Esso Osaka is a 280.000 dead-weight-ton VLCC with principal characteristics listed 


in Table 2.2. Abkowitz conducted a series of experiments outlined in [3] in the late 1970's 


Table 2.2: Esso Osaka Principle Characteristics 


LBP 1066.3 ft 
Beam 173.9 ft 
Draft 92.8 ft 

A, (rudder area) 1289.67 ft? 


AR (rudder aspect ratio) 1.538 
| A 314,410 lton 
LCG (aft of midship) 25 ft 





aud early 1950s to evaluate tlie EKF technique as a candidate for identification of ship 
livdrodyvuamue derivatives. The study produced favorable results aud simulations using the 
ideutified coefficients matched the full-scale data quite well. The successful results presented 
in [3] formed the basis for choosing the EKF method in this study. 

The matlenatical calculations involved in formulating the EKF are quite involved. The 
potential for error, therefore. is rather high. Thus. in order to ensure that no mathematical 
errors had occurred. it was necessary to test the filter against a dynamical model with 
kuown parameters. Successful operation of the filter in identifying the known parameter 
values would indicate that the filter had been properly formulated. The values of the 
coefficients appear in Table 2.3. The negative signs have been included in the coefficient 


definitions of Y, aud N, to follow the convention outlined in [49]. 
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Table 2.3: Esso Osaka Nondimensional Hydrodynamic Coefficients 


ma. — NS | 8.093e-4 






The results of the identification appear in Figures 2-3 and 2-4. These figures clearly 
indicate that the filter operates correctly. Notice how the parameter estimates converge to 
the exact values of the actual coefficients over time in Figure 2-4. This is indicative of the 
fact that the dynamical model in the filter is identical to the dynamical model producing 
the measurements. lhe tact that there are no modeling errors removes the requirement for 
process noise (.e. Q = 0). Thus. the filter trusts its own state estimates and converges 
quickly to the proper values based upon the initial measurements. Figure 2-3 shows how the 
estunates of the physical states converge immediately. This is due to the exact dynamical 
model in the filter and the fact that each of the four physical states are measured. Figure 2-4 
indicates that each of the parameter estimates requires about 3 minutes to converge. This is 
because each of these states must be estimated by the filter and some dynamic information is 
required before the filter can produce good estimates. After about three minutes, the error 
between the actual and estimated state vectors has decayed to zero. Thus. the computed 
salns are very simal and the filter has effectively “learned” the dynamics of the system and 


no longer requires measurement Information to produce good estimates. 
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Figure 2-3: State Estimates of the Esso Osaka During a 10°/10° Zig-Zag Maneuver 
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Figure 2-4: Parameter Estimates of the Esso Osaka During a 10°/10° Zig-Zag Maneuver 
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Results of simulations performed using the identified coefficients appear in Figures 2- 
5 aud 2-6. The identified coefficients reproduce the states and trajectory with very sinall 
error. This. again. can be attributed to the fact that the filter contained an exact dynamical 
model during the identification process. Note also that the trajectory shown in Figure 2-6 


in ideutical to the maneuver used to identify the coefficients. 
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Figure 2-5: Simulated States of the Esso Osaka During a 10°/10° Zig-Zag Maneuver 
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Figure 2-6: Simulated Trajectory of the Esso Osaka During a 10°/10° Zig-Zag Maneuver 
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Figure 2-7 shows the simulation results in a 10% steady turn maneuver. The simulation 
results are still quite accurate despite the differing maneuver. The error shown in the figure 
is less than 5% in the turning diameter. The error can most likely be attributed to roundoff 
errors in determining the coefficient values. This small error. however, is certainly ideal for 


control system design. 


Esso Osaka 10 Degree Steady Turn Maneuver 
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Figure 2-7: Simulated Trajectory of the Esso Osaka During a 10° Rudder Steady Turn 
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Chapter 3 


Nonlinear Simulation Model 


Many of the intricate details of the hullform described in the sequel have been onutted. This 
work uses a model developed in [43] to generate the data for the full-scale ship. Classification 
restrictions prevent. verification of this model against any full-scale data for the DDG-51. 
Thus. no attempt has beeu made here to do so. The focus of this work hes in using the 
SAEKF process to identify the unknown parameters in the equations of motion and develop 
control systems based upon the identified parameters. The goal then becomes using the 
identified values of the parameters to reproduce ship trajectories generated by the model 
from 43]. 

The sequel develops the general equations of rigid-body motion and the assumed form of 
the simulation model to be used in the SAEKF. Note that the form of the simulation equa- 
tions developed in this study are in no way related to the form of the equations used in (431. 
This ensures that the SAEKF does not merely estimate the parameters of a model identical 
to that which it coutaims. In fact. the form of the simulation equations is unimportant as 
long as the trajectory of the actual data is reproduced. Thus. the issue of coefficient can- 
ceHation described in [3] becomes moot. The goal of this study is not to identify individual 
coefficient values. but to reproduce trajectories. This imphes that the forces and moments 


have been adequately modeled. which is the most. important factor in control system design. 


3.1 The Arleigh Burke Class Destroyer 


The United States Navy (USN) Arleigh Burke (DDG-51) Class destroyer represents the 


state of the art in operational warships today. It 1s a twin-screw vessel powered by four 
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LAI2500 gas turbine engines coupled to controllable-reversible pitch (CRP) propellers. The 
vessel las a wider bean: than conventional warships, but is still considered to have fine 
lines. Table 3.1 lists the ship’s principal characteristics. The vessel uses two mechanically- 
linked rudders located downstream of the propellers to maintain course control. The ship's 
propellers maintain 100% pitch for all speeds above approximately 10 knots. and thus will 
be considered fixed for the duration of this study. The ship has port-starboard symmetry 


as do nearly all ships in existence today. 


Table 3.1: DDG-51 Principal Characteristics 


Length (LBP) 466 ft 
Beam (B) i 59 ft 
Draft (T) 

Displacement (A) 








Long Ctr Gravity (LCG) ; 2.8 ft aft midship 
Block Coeff (C,) | 0.922 
Prismatic Coeff (Cp) | 0.615 


3.2 Governing Equations of Motion 


3.2.1 Rigid-Body Inertial Forces and Moments 


The following sections develop the equations governing the inertial forces and moments ex- 
perienced by the body in dynamic motion. Simphfications ultnnately reduce the generalized 


equations to only those governing motions in the horizontal plane. 


3.2.1.1 Generalized Equations of Rigid-Body Motion 


The complete derivation of the generalized equations of motion for a rigid body nıoving on 
the surface of the carth begins by choosmg an inertial coordinate system with its origin 
located at the earth's center of gravity. This is necessary to account for centrifugal forces. 
coriolis forces, etc. The body experiences these forces due to the earth’s gravitational 
field and its relative hnear and angular velocity to that of the body. The effects of these 
phenomena. however. can be considered small (and thus neglected) for large objects such 


as ships. Reference [2] details this portion of the derivation, so it will not be reproduced 
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here. 


Based upou the assumptions outlined in [2]. the coordinate system is transformed into 


a body-fixed reference frame with 1ts origin located at the midship section. 


illustrates the coordinate system in a body-fixed reference frame. 


In the most general case. the 
ship may exhibit coupled mo- 
tious in all six degrees of free- 
dom shown in Figure 3-1. All 
three forces and all three mo- 
ments shown m Figure 3-1 act 
on the ship during normal op- 
eratiug conditions while under- 
wav at sea. Equating the forces 


and moments acting on the ship 


Motions: i 
u = surge | 
v= war Sa 
w = heave > 
p=roll 

q = pitch 
r= yin 






O = ortgin (midship) 


X = suree force 

d = rudder angle 
A $ TO 

E Z = heave force 

d N K = roll moment 

M = pitch moment 


N = yaw moment 


Figure 3-1: Body-Fixed Coordinate System 


Figure 3-1 


Farces and Moments: 


to the rates of change of hnear and angular momentum experienced by the ship yields the 


ecneralized equations of motion shown in equations 3.1 through 3.6. 


ea i +qw—rv— TO Irs IN PO (3.1) 
yo fi +ru — pw — yylr? +p°) + A a j) (32) 
Z =mi G + pv — qu — ao e a S nd +p) (33) 
Kpt grm fual w + pv — qu) — za(ù + ru — pu) (3.4) 
Af = I,q + (ls — an er, O + pv qn) (3.5) 
De, rpm al O II O - ro) (3.6) 


where 1 represents the mass of the body and Iz, ly, and 1, represent the moments of mertia 


about the appropriate axis. Because the longitudinal center of gravity 1s sufficiently close 


to the origin. and the ship is nearly symmetric, the cross-coupling inertia terms can be 


considered small and neglected. 


29 





3.2.1.2 Equations of Motion in the Horizontal Plane 


Despite the appearance of the tight coupling of all the motions im equations 3.1 through 
3.6. motions in the horizontal plane exhibit weak coupling to the out-of-plane motions [37]. 
49]. Therefore. motion in the horizontal plane can be considered separately from motions 


out-of-plane. This simplification implies the following: 


i — E= © on 


For a ship with no list (the usual case). the transverse component of the center of gravity. 
vo ınmst lie somewhere on the x-axis. Orthogonalitv of the coordinate svstem implies the 


following: 


Yo = 0 (3.8) 


Substituting equations 3.7 and 3.8 into equations 3.1, 3.2, and 3.6 yields the equations of 


motion 1 the horizontal plane: 


A =m i — rv — zgr?! (3.9) 

Y= n è + ru + 2 (3.10) 
r 

Ned rm E + ru) (3.11) 


The rıglt-haud sides of equations 3.9 through 3.11 describe the inertial forces aud mo- 
ments acting on the body during dynamic motion in the horizontal plane. The remaning 


derivations address the external forces that balance the inertial forces. 


3.2.2 Balancing Forces and Moments 


Many different forms for the equations governing the balancing forces and moments exist 


in the open literature. These models can be classified in three main categories [10]: 


p 


. Input-Output models. 


bo 


. Holistic (regression) models. 
3. Modular models. 
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Input-output models describe the direct effect of varying control parameters on the maneu- 
vering response of the ship. Holistic models treat the ship as a complete entity with forces 
aud moments described by a Taylor series expansion containing the pertinent kinematic and 
eeometric parameters. Modular models treat all major contributing elements as separate. 
interactive modules which can be developed and tested separately. Developments over the 
last decade indicate that the flexibility of the modular approach makes it very attractive 
for future applications [10]. Thus, this type of approach has been adopted in this work. 
The following sections develop the equations governing the forces and moments! which 
balance the inertial forces and moments derived in 3.2.1. The total force and moment can 


be subdivided into forces and moments resniting frem the following modular contributions: 
1. External influences (wind. waves. etc.), Xezt. Yert. Nezt- 
2. Steady-state effects. Xq, Yo, No. 
3. Propulsive devices (propellers. thrusters. etc.). Xp, Yp, Np. 
4. Rudder forces. Xp. Yp. Nr. 
5. Interactions between the fluid and the hull. Xp, Yp, Np. 


Thus. the total forces and moments may be described by the following equations: 


A = Nat = Xo ie Xp == XR + Xf (312) 
r Kot ne aa (3.13) 
N = Nen r Noti Ne aa N; | (3.14) 


3.2.2.1 External Forces and Moments 


This study considers maneuvering only in calm. deep water with no wind or current effects. 
This is not a likely operational environment for a full-scale ship. The goal of this study, 
however. is to reproduce maneuvering trajectories produced by tlie model developed in 
[43]. No environmental forces were considered in that study and therefore will not be 
considered here. Once the SAEKF method proves effective in the absence of environmental 
forces. the next logical step would be the addition of environmental forces to deterunue 


'X. Y. and N from equations 3.9 through 3.11. 
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their effects. Therefore. Xe. Yer. and Nex; are taken to be zero for the duration of 
the identification portion of this work. These effects will be considered, however, during 
simulation of controller designs. This is not a limitation because a scale model may be 
tested in a pool or tank where the environment can be carefully controlled to perform the 


identification. 


3.2.2.2 Steady-State Effects 


Steady-state effects consist of forces and moments present in steady-state motion. Thus, the 
steady-state Yo is taken to be the ship resistance at the steady-state forward speed when 
alh other dvnamne terms are zero Yo and No represent steady-state sway foree and vaw 
moment. respectivelv. These phenomena exist primarily on smgle-screw ships and mani- 
fest themselves through a tendency for the ship’s stern to “walk” in a particular direction 
when the propeller thrust is small. DDG-51 has two shafts and two propellers that rotate 
iun opposing directions. Thus. for tlis platform. the steady-state sway forces and yaw mio- 
ments cancel aud can be considered zero. Furthermore, because the steady-state forces and 


moments cancel. the dynamic forces and moments will cancel as well. Thus, 


Ya = E S (3.15) 


3.2.2.3 Propulsive Forces and Moments 


When a slip propels itself through water. the longitudimal force must be equal to the 
difference between the hull resistance and the propeller thrust. The open-water tlirust 


provided by the propeller can be expressed as follows [49]: 


To = pn DIT (3.16) 
(1 — w)u 

Jes AT 

nD om 


where n represents tlhe propeller speed. D represents the propeller diameter. and J represents 
the advance coefficient. The constants (Ko. Ky, and Ko) represent the coefficients in a 
parabolic fit of the thrust coefficient to the open-water propeller curve and {1—w) represents 


the Tavlor wake fraction. Inserting equation 3.17 mto equation 3.16 and carrying out the 
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requisite algebra results in the following expression: 


T, = Xp = mn? + mnu + nzu? (3-18) 
where: 
m =pD’Ko (3.19) 
m=(1-w)pD’K, (3.20) 
m = (1 - w)"pD*K, (3.21) 


Modeling of propeller thrust and torque introduces an additional state variable, n, into 
the equations of motion. Appropriate models for this variable depend heavily upon the 
type of propulsion imachinery the plant contams. The DDG-51 propulsion system consists 
of four LA{2500 gas turbine engines mechanically coupled to two shafts via a set of reduction 
gears. Reference [49] reconmends the following form for modeling a gas turbine propulsion 


system: 


a nrAngQe — QP 


50) 
n 271, Sg 


where Qe and Q, represent the engine and propeller torque, respectively, [, represents 
the polar moment of inertia for the entire engine/gear/shaft/propeller arrangement. À rep- 
resents the reduction gear ratio. and 7, and ng represent the propeller relative rotative 
efficiency and reduction gear transmission efficiency. respectively. This model is intended 
ouly to adequately model the changes in propeller speed during maneuvers. It is not in- 
tended to completely capture the performance of the engine itself. Adequate modeling of 
propeller speed changes captures the salient dynamics required for maneuvering and tra- 
jectory simulation. Therefore. this model suffices for the task at hand. Thus, the engine 


model used in this study is the same as that used in [43] °. 





QE = Qmar | (afr +b) = cf (320) 
Oe = pn’D’Ko(J) = pn? D° lQ. +QıJ + 027 (3.24) 


“Reference [1] provides a more detailed model to capture the performance of the engine. 
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where all coefficients in equation 3.23 are taken directly from [43]. f, represents the fuel 


rate. and the constants (Qo, Q1. and Q2) represent the coefficients in a parabolic fit of the 


torque coefficient to the open-water propeller curve. 


3.2.2.4 Rudder Forces and Moments 


Modeling the rudder forces and moments is somewhat more complicated thau the derivations 


presented thus far because the angle of attack decreases due to the sway and yaw motions. 


Figure 3-2 illustrates the coordinate system and vector diagrams of the rudder forces. The 


ultimate goal is to determine expressions for Xr and Yr. An expression for Nr then follows 





Ap= Rudder Area 
Ci (X R) = Lift Coefficient 
CA R) = Drag Coefficient 





Figure 3-2: Vector Diagram of Rudder Forces 


It is well-known that rudders are foils that act as lifting devices when positioned at some 


augle of attack. ar. measured with respect to the oncoming flow velocity. Ur. The finite 
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RH Ko 


radius at the leading edge of the foil requires that the lift force, L. act. perpendicular to 
the oucoming flow [49]. Orthogonality of lift and drag requires that the drag force. D. act 


parallel to the flow. The lift and drag forces are defined in equations 3.25 and 3.26. 


1 

L= 5PARURCL(aR) (3.25) 
l 

D = 5PARURCD(aR) (3.26) 


When v aud r equal zero. œp is simply the rudder angle, ô. When the ship experiences 
sway and yaw iotions. however. ag decreases as shown in Figure 3-2 3. The angle w 
represents the decrease in the augle of attack due to swav and vaw. Thus, the instantaneous 
hft and drag forces depend upon the instantaneous flow velocity. Up. and instantaneous 
angle of attack. ar. The quantities shown in Figure 3-2 are determined froin the relations 


shown in equations 3.27 through 3.29. 


Le =u + (v+ rr} (32T 

ar =ð +4 (3.28) 

p= arcan (=) (3.29) 
u 


The lift coefficient. Cr (ar), and drag coefficient. C'p(ar). can be approximated with any 
appropriate model. A typical value for the drag coefficient of a foil is 0.0085 [31]. The Hift 
coefficient can be modeled as a constant for small angles of attack [49]. Use of the modular 
rudder model allows tle use of tlis approximation even ni maneuvers where it appears to 
be invalid (1e. a 20° rudder turn). This is because the actual angle of attack decreases 
rapidly witl: respect to the nominal as sway speed and yaw rate increase. Reference [3] used 
a 20° rudder turn to validate nonlinear coefficients identified in more violent maneuvers. 
This can be taken to mean that in that work, this type of maneuver was considered to be 
at least nuldly nonhuear. The sequel will show. however. that this approximation works 
quite well for maneuvers of this type. For foils with aspect ratio (AR) greater than 1, [49] 


“Figure 3-2 shows an apparent increase in the angle of attack for illustration purposes onlv. « always 
opposes 6 to allow the forces to balance in steady state. 
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provides the following expression for the linear lift coefficient, Cz: 


Cr = St = ( Me (3.30) 
“dar 0.187 TAR MAR)? | 

Thus. taking Cp = 0.0085 and computing Cz, from equation 3.30. the lift and drag forces 
may be computed from equations 3.25 through 3.29. Once the lift and drag forces are 
known. the rudder forces are coniputed with equations 3.31 through 3.33 where the factor 


of two accounts for the fact that the ship has two rudders. 


Xp =2(Lsiny — Dcos Y) 3) 
Yr = —2(Lcos y + Dsmyp) (3.32) 
Nr = Tp (3.50) 


3.2.2.5 Hull/Fluid Interaction Forces and Moments 


The forces and moments generated by the interaction between the fluid and the hull are 
functions of several variables. Thus, Xp, Yp. and N; may be described by the following 


equations: 


Af = I a, T} (3.34) 
a= joven) (3-35) 
Ny = falu.v.r) (3.36) 


The open literature provides many expressions for f1. f2, and f3 along with detailed deriva- 
tious and underlying assumptions?. The derivations and assumptions will not be repeated 
hiere with the understanding that all equations presented in the sequel conform to all as- 
sumptions described iu the references from which they were taken. 

The major differences between the proposed models he in their treatment of the non- 
lincaritics. Abkowitz proposed the most well-known form of the equations outlined in [37]. 
This model proposes a third-order Taylor series expansion of the fluid forces and moments 
about some steady-state forward equilibrium speed. Hwang [24] further refined this model 
based upon full-scale maneuvering trials conducted on a series of large tanker ships [3]. 


ISee 223. 7). [131. 28 [52}eand [54]. 
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Blauke 5. proposes a simpler form of the equations intended to capture the most important 
nonlinearities in terms of speed and propulsion loss. The form of the simulation equations 
chosen for this study is a hybrid of the models proposed in [5] and [24]. It should be noted 
that each of the models mentioned here are based upon the holistic method described in 
3.2.2. Because this work employs the modular method, the hydrodynamic derivatives asso- 


ciated with the rudder forces in these models? are captured by the rudder model described 


Ap = filucv.r) = Xyu+Xpt+ Xp (3.31) 
Dar yO Yr HY, (3.38) 
Nr = falu,v,r) = Nov + Nar + Ny, (3.39) 


Added mass and added inertia terms associated with accelerations (i.e. ù, ©, etc.) 
can be calculated to within sufficient accuracy for purposes of control systems through 
hvdrodvnamic strip theory. The terms Xn. Yp. and Np in equations 3.37 through 3.39. 
however, represent some combination of unknown damping and nonlinear terins whose 
values no existing hydrodynamic theory can predict witlr any level of accuracy. These 
terms are currently determined through extensive model testing. Thus, the linear damping 
and pertinent nonlinear terms must be estimated to provide at least some knowledge to be 
beneficial in the design of a control system. 


3.3 Nonlinear Simulation Equations of Motion 


Completing the steps that follow results in the final form for the nonlinear simulation 


equations of motion. 
l. Apply the simplifications and assumptions discussed in 3.2.2.1 through 3.2.2.3. 
2. Introduce the resulting expressions for X, Y, and N to equations 3.9, 3.10, and 3.11. 


3. Solve the resulting system of equations for the state derivatives and obtain a system 
i the state-space form of equation 2.10. 


Pe eto: 
“Equations 3.37 through 3.39 employ standard SNAME notation for hydrodynamic derivatives. Thus. 
the coefficients in the Tavlor series expansion =: L, etc. are absorbed by the hydrodynamic coefħcients. 


See [37] for details. 
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The final model which results froın completing each of the steps listed above appears ın 
Table 3.2. Coefficients of added mass and added inertia terms in Table 3.2 are determined 
from valid hydrodynamic theory and knowledge of the ship geometry. Coefficients of the 
linear damping terms and any nonlinear terms (i.e. components of Xa. Yn. and Np) are to 


be treated as the unknown parameters to be estimated using the SAEKF. 


Table 3.2: Nonlinear Simulation Equations 





| U — = 

| m-X; 

| , — Uz — Ni) fo — (may — Yi) fs 
| 2 

| ha 

| . _ NrAQE - Qp 

| | er, 

| 

| 


= nn? + nnu + nzu? + X,+ Xp 


Js = Np +N R 


p 


pet, — N; (n Yo) m n a N) 


f= Vn + Yr | 
| 
| 
| 
m=(1-w)pD’Kı | 


SS pD* Ko 
| 
| 
| 
| 


13 = (1 — w)? pD? Ko 


Coefficients of inertial terms such as };; and N; to be determined from hydrodynamic theory. 
Pertinent terms in X,. Y), and N, to be estimated with SAEKF. 
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Chapter 4 


DDG-51 System Identification 


Effective control system design relies upon simplified models that adequately describe the 
dominant dynamics of the plant to be controlled. These simplified models inust contain the 
minimum number of parameters required to sufficiently describe the forces and moments 
acting on the plant. The equations developed in Chapter 3 are intended to describe thie 
complete noulmear plant dynamics. The dominant forces and moments are often captured 
by linear terins over a surprisingly wide range, however. Modern control system theory 
provides the capability to design controllers that are robust in the face of modeling errors. 
Thus. the first step in design is to design a controller based upon a linear model and test 
its performalce on the nonlmear plant. 

The linear ship dynamics are governed by added mass and added inertia terms. as 
well as linear damping terms [37]. Hydrodynamic theory provides sufficiently accurate 
methods for calculating the former. Plant stability, however, 1s governed by the latter, for 
which no theory provides sufficiently accurate results. These terms must be determined 
experimentally. Several experimental methods for determining these terms are outlined in 
137). The drawback to these experimental methods is that they require expensive equipment 
and labor in addition to a scale model to collect the necessary data. A method of determining 


these terms using only the scale model is one of the focal points of this study. 


4.1 Linear Ship Dynamic Equations 


Ships operating at sea most. often conduct maneuvers that le within the linear regime. For 


example. slips in normal operating conditions do not usually use very large rudder angles 
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m turns. The most common rudder angles are often between 10-15°. This type of maneuver 
lies well within the lincar regime. During such maneuvers, forward speed loss and propeller 
rotational speed loss are small. They do exist, but for the purposes of control system design. 
they are sufficiently small such that they do not significantly effect the performance of the 
controller. This implies that. the surge equation, ù, and the propeller speed equation, n. 
in Table 3.2 are uncoupled from the sway equation. ù, and the yaw equation, r. Thus, the 
identification process reduces to the determination of fo and fs. 

Table 3.2 shows that fo and fs each contain two terms which describe the hull damping 
forces and the rudder forces. The rudder force model developed in 3.2.2.4 is based solely 
upon wine theory and the kinematics of the problem. Thus. this force is assumed to be 
accurately modeled aud requires no identification. Therefore, the identification problem 
rednces to determining the hull damping forces, Yp and N). Abkowitz [37] suggests that 
m the nonlinear case. these terms are adequately modeled with a third order Taylor se- 
ries expansion. This expansion. however, contains a very large number of terms which. is 
intractable from a control system perspective. For linear maneuvers, however. the expan- 
sion may be truncated after the first order terms. This truncation leads to the following 
expressions for fọ and fs: 


fa = -Y,v - (m - Y.)r a) 


y 


f3 = —N,v — [mzxg — N,)r (4.2) 


Thus. for the case of linear maneuvers, the identification problem is reduced to the deter- 
mination of the parameters —Y,, (m — Y,), —N,. and (mxg— N,). The negative sign on Y, 
aud N, was introduced to maintain the convention outlined in [49] that —Y, should have a 


large positive value. 


4.2 Initial Identification 


The 10°/10° zig-zag maneuver proposed in [3] was used to determine the linear coefficients. 
This type of maneuver lies well within the linear regime and also provides “persistence 
of excitation”. Reference [34] states that “An open-loop experiment is informative if the 


input is persistently exciting.” This effectively states that the system dynamics must be 
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continually excited in order to get any information about the parameters. A steady turn. 
for instance. would quickly reach a steady state and provide no more dynamic information 
to update the filter. Thus. persistent excitation is crucial to successful identification. 

Recall that the model generating the measurement data was taken from [43]. This nıodel 
was adapted from a form proposed by Inoue in [26]. The model is parametrically based on 
full-scale data obtained from experiments performed over a range of operating conditions. 
The model suffers. however, during a simulated zig-zag mancuver. The implementation of 
the model in |43] suffers from a drift effect due to the rudder model. Figure 4-1 illustrates 
this effect. 
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Figure 4-1: 10°/10° Zig-Zag Maneuver Using the Rudder Model in [43] 


Initial attempts at identifying the coefficients using this model produced biased results. 
Because no a priori information was available for the value of the unknown parameters, the 
ouly way to determine success or failure of identification was through simulation. Figures 


4-2 aud 4-3 show the values obtained through identification and the resulting simulation. 


4.2.1 Bias and Divergence in the Extended Kalman Filter 


Divergence and bias in the Kahnau Filter is well-documented in the open literature. Diver- 


vence aud blas. however. do not imply filter instability. The EKF possesses some guaranteed 
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Figure 4-2: Parameter Estimates Using the Rudder Model from [43] 
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Figure 4-3: Trajectory Simulation Using the Rudder Model from [43] 
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stability properties under a few mild conditions [20]. Thus. the filter can exhibit vappar- 
ent divergence” [12] and give incorrect (biased) parameter estimates. References [12] and 
44] propose that bias is caused primarily by errors in the filter model. This suggestion is 
well-supported by [14] and [30]. Reference [22] provides numerical examples of the effects 
erroneous models can have on the performance of the filter. Each of these suggest that 
the model of the plant used in the filter must be sufficiently close to the actual plant to 
euarantee rehable results. 

The process noise covariance matrix. Q. used in the filter model is meant to account for 
modeling errors. Thus. the choice of Q represents a critical filter design parameter. If tlhe 
process noise ts chosen too small. the filter “learns” the wrong state values too well. This 
means that the error covariance decays too rapidly and the filter ignores any additional 
information contained in the measurements. Thus. the process noise may be said to “drive” 
the filter aud ensure that it does not place undue weighting on its own estimates. Increasing 
the process noise. however. lias a strong effect. on the convergence of the filter. Thus, it can 
not be chosen too Ingh or the filter will not converge. In an effort to account for the bias. 
identifications were attempted using a wide range of process noise values ranging from 0 
aud covering two orders of magnitude. Similar results were obtained m all cases. 

Much work has gone into developing an analytical method for deternnning the process 
noise covariance. Refereuce [46] proposed a method for handling noise covariances through 
post-processing of data. References [18]. [39]. and [40] all propose methods for estimation 
of the covariances in real-time. Friedland proposes another method for estimating the bias 


terms in real-time in [14]. 


4.2.2 Treatment of the Biased Estimates 


Favorable application results were demonstrated in each of the works cited above. Ap- 
plication of these techniques to the ship maneuvering problem, however, has not been 
demonstrated to the author’s knowledge. The most favorable application of the SAEKF 
technique to the ship maneuvering problem is summarized in [3], which is the culmination 
of 10 vears of research conducted at the Massachusetts Institute of Technology (MIT). The 
details are outlined in the combined works of [6]. [21], [24], [35], and [47]. It is evident in [24] 
that bias was a problem as well. although it is not explicitly described. In that work, the 


assumed form of the simulation model was significantly altered from the Taylor expansion 
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form: proposed in [2] and (37] to obtain favorable results. 

This rationale has been adopted in this work as well. The drift illustrated ın Figure 4-3 
most likely results front numerical problems in the model from [43] due to its parametric 
nature. The linear dynamics of ship maneuvering are well-documented and the assumed 
form shown in equations 4.1 and 4.2 are not likely to be overly erroneous. For these reasons, 


the rudder model used in [43] was replaced with the rudder model derived in 3.2.2.4. Any 


further references to this model assume this modification has been made. 


4.3. Identification of the Linear Damping Coefficients 


4.3.1 Noise Parameters 


The magnitude of the process noise covariance matrix, Q. and the sensor noise covariance 
matrix. R. determine the weighting the filter applies to its own estimates and the measure- 
ments. respectively. During the course of the study, the trade-off between the two became 
evident. Given the sophistication of current digital modern measurement devices. the non- 
dimensional measurement noise was chosen in all cases as R = 0.011}. When the process 
noise is too small. the filter“trusts” its own estimates too much and, thus. gives a very low 
weighting to the measurements (i.e. low gains). This introduces the bias discussed in 4.2.1. 
Figures 4-4 and 4-5 illustrate this phenomenon. 

The identification uses the nonlinear model from [43] to generate the data and the filter 
emplovs the linear nodel described in 4.1. Figure 4-4 shows that the states are apparently 
tracked quite closely. Note. however. that the filter never updates the surge speed. u. This 
is caused by the fact that the absence of process noise implies perfect modeling of the plant 
dynamics. Therefore. the filter has no reason to update its own propagated estimates (1.e. 
low gains). The close tracking of the sway speed, v, and the yaw rate, r, indicate that the 
uonlinear terms in these equations contribute very little during a linear maneuver. This 
indicates that the dynamics are modeled with sufficient accuracy. Thus, the filter is justified 
im trusting its own propagated estimates. 

Figure 4-5 shows the parameter estimates during the identification. Note how each of 
the parameters converges nicely to a constant value. This indicates that the error covariance 


'This value is of the same order of magnitude as the non-dimensional forces, which is believed to be 
conservative. 
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State Estimates with Q=0 
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Figure 4-4: State Estimates During a 10°/10° Zig-Zag Maneuver with Q=0 


Parameter Estimates with Q=0 
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Figure 4-5: Parameter Estimates During a 10°/10° Zig-Zag Maneuver with Q=0 
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has decaved to zero aud the filter stops updating the estimates. Recall the expression for 


the ganis. Ky. from Table 2.1. This equation imphes the following: 


Ky, = 0.as Py(—) = 0 = 3, (+) = &(—) (4.3) 


The results in Figure 4-5 are inisleading, however. since they exhibit apparent, divergence 
112]. Figure 4-6 illustrates the results of the simulation after identification. The steady turn 
is chosen as the simulation maneuver to illustrate the performance of the filter estimates 111 


a maneuver other than that used to identify the coefficients. Figure 4-6 shows that. while 


Simulation After Identification with Q=0 
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Figure 4-6: Simulation of a 10° Rudder Steady Turn with Q=0 


the filter does converge. the estimates produce an unstable plant. 

If perfect modeliug of the system dynamics were possible, and the estimates still proved 
to be biased. the designer may introduce some “fake” noise into the system. This “fake” 
noise would continue to drive the parameters once the state estimates have converged and 
cause the filter to weight the measurements more heavily. This technique is outlined in [18] 
aud suggests that a value of 0.1% of the nominal parameter value is suitable for this purpose. 
This approach fails. however, for two reasons. First, using a “nominal” parameter value 
unplies some a priori knowledge of the plant. In this case, as in most cases in the physical 


world. this does not exist. Second. because the non-dimensional parameter values are small, 
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once the state estimates have converged the weighting on the measurements will be very 
sulall as well. Thus. the small weighting on the measurements tends to drive the filter very 
slowly. Thus. it will not converge in any reasonable period of time. Figure 4-7 illustrates 
the effect of the “fake” noise. Note how the parameters appear to be slowly decaying with 
time. The slope is so gradual, however, that the filter will require an unacceptable amount 


of time to converge (if it does at all). 


Parameter Estimates with Q=0 and "Fake" Noise 
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Figure 4-7: Parameter Estimates 


Given the previous considerations. the appropriate noise parameters must consist of 
some combination of process noise and “fake” noise such that the filter continues to utilize 
the information contained m the measurements. The appropriate combination of these noise 
parameters was chosen through an iterative process of identification and simulation. Large 
process noise was required to effectively track the states. This large process noise, however. 
caused rapid convergence to the true state values. Thus, “fake” noise was also required to 
continue to drive the parameters to their true values. The combination of the two noise 
parameters allowed the filter to update its dynamical model, as well as obtain information 
from the measurements until a balance was achieved. The measurement noise was chosen 


as described in 4.3.1. The final values of the noise parameters are listed in Table 4.1. 
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well. The filter also does an excellent job of tracking the sway speed. v. and vaw rate. r. 
This. again. is due to the fact that sway and yaw are measured quantities and the process 
noise is high. Further, because the nonlinear effects are small in this type of maneuver, tlie 


dynamic model for tliese two ternis Is quite accurate. 


4.3.3 Parameter Identification 


The four additional states in the filter represent the unknown parameters to be estimated. 
Previous sections demonstrated the filter's ability to track the physical states closely, wlule 
viving inaccurate (biased) parameter estimates due to modeling errors and erroneous as- 
Mens i tha noise statisties. The iterative process of identification and smimlation 


produced the parameter estimates illustrated in Figure 4-9. 


DDG-51 Estimated Hydrodynamic Coefficients During a 10°/10° Zig~Zag Maneuver 
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Figure 4-9: Paraineter Estimates During a 10°/10° Zig-Zag Maneuver 


With the exception of Y,. the parameter estimates exhibit variatious about some value. 
This could prove to be problematic if tle goal of the identification were to establish exact 
coefficient values. The goal of control system design is to develop controllers that perform 
well with only an approximate knowledge of the coefficients. In fact, parameters that niodel 
the state values within 40% of the actual values often leads to acceptable control system 
performance. Fortunately. tlis goal is often achieved. Thus, to determine the coefficient 


value to be used in the simulation. the mean value of the last 67% of the parameter estimates 
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was taken to be tlie parameter value. 


4.3.4 DDG-51 Simulation After Identification 


The ultimate goal of an automatic maneuvering control system on a slıip is to provide a 
rehable method of maintaining a desired trajectory. Linear control system theory has been 
proven to be quite reliable in designing control systems for plants in which nonlinearities are 
weak. Control of an inverted pendulum represents an excellent example of this assertion 7. 
Because ship maneuvers normally lie within the linear regime, where nonlinear effects are 
small, a control system design based upon a linear ship model could very likely perform quite 
wol’ Fine 10 and 1-11 illustrate the sinmlatión results using the ideutified parameto 
values. 


Simuiated States After Identification 
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Figure 4-10: Sunulated States During a 10° Rudder Steady Turn 


Figure 4-10 shows that the identified model tracks the actual nonlinear model quite 
closely iu sway aud yaw. It does not. however. track surge and propeller revolutions at all. 
This is because the simulation model iguores forward speed loss aud propeller speed loss 
(ie. ù = n = 0). This is a good assumption because the actual speed loss in the simulation 
is about 20% and propeller speed loss is approximately 10%. These represent ideal errors 
for the applicatiou of linear theory to a nonlinear physical system. The lack of speed loss 


“See example 2.1-1 in [33]. 
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Figure 4-11: Simulated Trajectory During a 10° Rudder Steady Turn 


accounts for the Ingler values of sway and yaw. Higher speed in a turn imphes a tighter 
turn. and thus. higher values of sway and yaw. 

Figure 4-11 shows the trajectory produced by the identified model. The error in the 
turning diameter between the identified and actual models is approximately 15%. This 
situation is ideal for the application of linear control theory to develop automatic controllers 
for the ship. Note that. the identified model has a tighter turning diameter than the actual 
model. This was expected due to the higher values of sway and yaw shown in Figure 4-10. 


The identified model maintains higher speed throughout the turn. Thus, it has a smaller 


turmime diameter. 


ol 








Chapter 5 


Controller Design 


The previous chapters focused on the identification of the linear damping ternis in the 
DDG-31 equations of motion. The identification process resulted in determination of the 
hydrodynamic coefficients that produce the dominant forces and moments during hnear 
ship maneuvers. Thus, these terms may be employed to develop simple models which may 
form the basis for simple course-keeping controllers. 

The sequel deseribes the design of two controllers based upon simplified ship models. 
The first controller is a heading autopilot. The second controller regulates the lateral 
deviation (cross-track error) from the desired reference trajectory. The desired frequency 
responses are developed using a loopshaping technique. The desired control laws are then 
determined from the design loopshape. Simulations are then performed on the nonlinear 


ship model using the hnear controllers to evaluate their performance. 


5.1 Controller Design via Loopshaping 


Controller design via loopshaping is well-documented in many texts on feedback and optimal 
coutrol system design. More detailed versions of the following presentation are described in 
111]. [33]. and [53]. 

The overall goal of controller design via loopshaping is illustrated in Figure 5-1. The 
Loop Gain should be designed such that the closed-loop plant exhibits robust performance 
111]. This means that tle entire system (plant and controller) exhibit good performance 
and tracking at low frequencies within the system bandwidth. The system must. also be 


robust to high-frequency disturbances such as sensor noise and unmodeled dynamics. The 
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Figure 5-1: The Goal of Loopshaping 


sequel describes how this process 1s accomplished by the control system engineer. 


5.1.1 Feedback Control 


Figure 9-2 represents a typical feedback control system with the plant disturbance input. 
d. reflected to the plant output as in [33]. Reference [11] considers the feedback loop with 
the disturbance reflected at the plant input. The derivation of [33]. however. makes the 
effect of the disturbance and desired properties of the loopshape a bit. more clear. Thus. it 


is primarily that derivation that will be presented here. 





Figure 5-2: Typical Feedback Loop 


The goal of a feedback control system is for the plant output, y, to asymptotically track 
a destred reference output, r. The plant, P, however, is acted upon by disturbances, d. 
which represent external disturbances and modeling errors. Further, the measurements of 


the output are influenced by sensor noise, n. The goal of the controller, C, is to compute 


93 





the control input. «. that regulates the output and maintains the desired reference value 


despite the disturbances and measurement noise. Defining the tracking error. e. as follows. 


e=r—y (aal) 


exammation of Figure 5-2 gives rise to the following expressions: 


Very on (5.2) 


y = PCu=d (5.3) 


9.1.2. System Sensitivity, Cosensitivity, and Loop Gain 


To illustrate the ideas behind loopshaping theory. one requires expressions for the closed- 
loop transfer function relations of the desired outputs. y and e. to the deterministic and 
stochastic inputs. r. d. and n respectively. Substituting equation 5.2 into equation 5.3 


results m the following expression for y 


Y => PONSA (5.4) 
(I+ PC)y= PC(r -n)+d (505) 
(A. HO ee Gea ee Glad (5.6) 
PONEN AAA (5.7) 


where the final operation is justified because PC is sanare and invertible [33]. The expression 


for the error. e. may be derived in a similar manner. 


e=r-PCv-d (5.8) 
e=r-PC(r-y-n)-d (5.9) 
e=r-d-PCie-n) (5-10) 
(I + PC)e = r — d + PCn (5.11) 
e=(I+PC)"(r-d)+(I+ PC) !PCn (5.12) 
e = (ME PCO dd) Beim PC) In allay 


Equations 5.14. 5.15, and 5.16 define the system sensitivity, S, cosensitivity, T, and loop 
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gain. L. respectively. 


Seen (5.14) 
TE POE BOY” (5.15) 
Die (5.16) 


Notice that S+7T = I. Thus. S and T may not both be made arbitrarily large. Substituting 
equations 5.14 and 5.15 into equations 5.7 and 5.13 results in the following expressions for 


the output. y. and the tracking error, e: 


y=T(r—n)+ Sd LOA 


e=S(r—d)+Tn (5.18) 


System disturbances are low-frequency phenomena such as wind and wave forces exerted 
on the ship hull. Examples of high-frequency phenomena are sensor noise and structural 
vibrational modes. The goal of the controller is to respond to low-frequency disturbances 
within the bandwidth of the system aid attenuate disturbances outside the system band- 
width where its behavior is not as well-known. For example. the controller should respond to 
wave disturbances to maintain the desired output. but should not respond to lugh-frequencv 
structural vibrational modes whose dynamics were neglected in the system model. 


At low frequencies. equations 5.17 and 5.18 may be expressed as follows: 


yx Tr+Sd (5.19) 


ex S(r —d) (5.20) 


At high frequencies the same equations may be represented as follows: 


because the refereuce signal is typically a constant or at least slowly-varying. Thus, the 
system sensitivity. S. should be small at low frequencies where the disturbances, d, are sig- 


nificant to ensure small tracking errors. Additionally. if S is small, T = 1 at low frequencies 





(because S + T = 1) such that y = r. Converselv. the system. cosensitivity. T. should be 
small at high frequencies where sensor noise and other high-frequency phenomena are sig- 
uificant. If these conditions are satisfied, the controller will exhibit robust performance [11]. 
This implics that the controller provides good tracking of reference inputs (performance) 
while sunultaneously rejecting high-frequency signals where the plant model is not so well 


understood (robustness). 


5.1.3 Performance Specifications 


Frequenev-dependent weighting functions are required to establish the low and high fre- 
queney performance bounds shown in Figure 5-1. These funetions establish the minimumn 
aud maximiun open-loop system gaims over all frequencies. Thus, one weighting function 
should be large (i.c. large positive gains) at low frequencies to ensure good performance. 
The other weighting function should be small (i.e. negative gains) at high frequencies to 
ensure robustness. 

Determination of these weighting functions. coupled with design of the loop gain. L, 
then allows development of frequency-domain performance specifications. The idea behind 


performance specifications may be illustrated on a Nyquist. plot [11] as in Figure 5-3. For 
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Figure 5-3: Nyquist Plot of Performance Specifications 





evod tracking. the loci of L must lie outside the disk of radius W’1 over all frequencies. To 
eusure robustness m the face of multiplicative plant uncertainties, the disk of radius W2 
must uot pass through the point —1 + 07 over all frequencies. 

The goal of the controller is robust performance which implies good tracking and distur- 
bance rejection at low frequencies and attenuation of high-frequency signals such as sensor 
nolse aud modeling errors. Recall that the system sensitivity function, S, should be small 
at low frequencies to provide good performance and the system cosensitivity function. T. 
should be small at high frequencies to provide robustness. Considering Figure 5-3 and the 
definitions of the quantities derived previously, it can be shown [11] that the frequency- 
domain performance specification that ensures robust performance is viven bv the following 


eCXpression: 
is) + peri] EE (5.23) 
DS 


5.1.4 Design Criteria 


The length between perpendiculars (LBP) for the DDG-51 class ships is 466 feet. The main 
low-frequency phenomena in the open sea that influence the ship motions consist of waves 
and current. Currents are most often constant or at least very slowly varying. Thus. waves 
represent the phenomena that regularly contain a frequency component. 

The frequency of open-ocean gravity waves is related to their characteristic wavelength. 
The motion of the ship is also related to this wavelength. According to [49], the closed- 
loop system should respond to wavelengths of the same order as the ship LBP. Thus. the 
closed-loop system should respond to all wave disturbances of length 400 feet or greater. 
This corresponds to a characteristic circular frequency of 0.7 radians/second (Period = 9 
seconds). To allow for some flexibility in the loopshaping, the design crossover frequency. 
w,. should he in the range between 0.5 - 1.5 radians/sec. 

Currents tend to “blow” the slip off course in the same way that wind blows an airplane 
off its intended course. If this disturbance is not accounted for, the ship would move far from 
its intended track. Thus. the ship must adjust its heading angle such that it accounts for the 
drift effect due to the current. Because currents are very low frequency disturbances. the 
closed-loop system must track them very accurately. Thus, the closed-loop system should 


track disturbances with an error, @. less than 1% at frequencies below 0.01 radians/sec. 


of 





This corresponds to open-loop system gains of at least 40 dB below 0.01 radians/sec. Thus. 
the open-loop system gains should lie above the iow-frequency performance bounds and 
below the high-frequency performance bounds illustrated in Figure 5-1 while achieving a 
reasonable closed-loop step response. The ultimate test of any controller design is how it 
operates on the actual rather than the nominal plant. The reasous for not ae time 


domain performance criteria will be addressed in the sequel. 


5.2 Autopilot Design 


The goal of a ship autopilot is to maintain the ship’s actual heading around some slowly- 
varying reference heading. The autopilot must perform tlis function m the presence of 
disturbances, sensor noise. and modeling errors described in 5.1. The design speed for all 
siinulations is 22 knots. This corresponds to the maximum ship speed in a plant config- 
uration with two of the four main engines on-line. one driving each shaft (split plant). 
Non-dimensionalization of the equations, however. allows this speed to be adjusted without 
affecting the results. The split plant configuration is required, however, to validate tlıe 
assumption in Chapter 3 that the steady-state forces and moments due to the propeller 


cancel. 


5.2.1 Dynamic Model 
Low frequency ship motions may be described by a first-order model attributed to Nomoto 
113]. This model may be represented as follows: 


Ty +4 = K8 (5.24) 


where y represents the ship heading angle and ô represents the rudder angle. Note that 
the sway dynamics have been removed from the equation of motion. T and K represent 
the ship time constant aud gain. respectively. These terms may be determined from tle 


hydrodynamic coefficients identified 1 Chapter 4 using the procedures outlined in section 


08 


5.3.2 of 13). This leads to the following values for T and K: 


T = 1.9287 (5.25) 


K = 2.04 (5.26) 


Thus. the plant transfer function. Py from input, 6, to output, w. appears in equation 5.27 


and represents the change m heading angle to a change in the rudder angle. 


Y K 


P = — SO 
"ai s(T's +1) 


(5.27) 


The ship's rudder contains its own dynamics. For instance, the rudder does not in- 
stantaneously change its angle when conımanded by the helmsman. According to [24]. the 
rudder dynamics do not play a major role in the identification process because they are 
much faster than the ship dynamics. For control purposes, however, the rudder dynanı- 
les do play a role. Classification restrictions prevent the disclosure of the actual DDG-51 
rudder swing rate. Thus. it was assumed in this work that the rudder dynamics could be 
modeled by the following transfer function. Psg, that maps the commanded rudder angle to 


the actual rudder angle: 





Thus. the new plant, P. may be represented as follows: 


P = Pa Pog (5.29) 


This new plant. P. corresponds to the plant block in Figure 5-2. To guard against. unrea- 
sonably large rudder angle commands during maneuvers, a rudder command angle limiter 
was added to the model as well. This limits the maximum allowable commanded rudder 
angle during maneuvers to 30°. This guards against the actual rudder angle approaching 
the physical stops on the actual ship. The goal now is to design the controler, C, using 


loopshaping to achieve robust. performance. 


99 





5.2.2 Loopshaping Controller Design 


[In some instances. it is possible to shape the loop gain. L. to obtain the desired frequency 


response. The controller may then be determined as follows: 


C= (5.30) 


L 
P 
In other cases. the controller must be designed directly to obtain the desired response. The 
autopilot design followed the former procedure. The open-loop plant, P, is straight-line 
stable [37] as are many ships. Thus, the controller is not required to stabilize the plant. 
To sive this a physical representation. imagine the ship without any contro] encounters 
a current. The current will cause an initial yaw rate, but the large wetted surface area in 
the stern will develop forces and moments that counteract the initial rate of yaw. Thus. the 
vaw rate will decay back to zero after a finite period. Because the heading angle is simply 


integral of the yaw rate. or 
br (5.31) 


the slap will turn to a heading that faces the current and stay on that. heading when r = 0. 
Thus. the goal of the autopilot must be to ensure that the slip continues on its desired 
heading despite the current. 

Thus. the loop gain for the autopilot design may be modeled simply as a second-order 


system with an integrator whose transfer function may be represented as follows: 


pe 
s(s* + 2¢w, + w2) eS 

wn = V60 (5.33) 
¢=2 (5.34) 


Figure 5-4 illustrates the results of the loopshaping design. The loopshape has gains of 
at least 40 dB at. frequencies below 0.01 radians/sec. It lies above the low-frequency per- 
formance bound indicated by W1. This indicates good tracking and disturbance rejection 
at low frequencies. The loop gains roll off at high frequencies indicating good attenua- 


tion of high-frequency signals and unmodeled fast dynamics. They further lie below the 
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Figure 5-4: Loopshape for Autopilot Design 


hıgh-frequeney performance bound indicated by Wa. The design also satisfies the robust 


performance requirement given by equation 5.23. 


An indication of overall plant stability may be obtained from the Bode plot. Figure 5-5 


indicates a gain margin of approximately 30 dB and a phase margin of approximately 64”. 


Reference [49] indicates that a good rule of thumb for gain and phase margius are 3 dB and 


30°. respectively. Thus. the autopilot design falls well within the acceptable range. 
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Figure 5-5: Gain and Phase Margins for Autopilot Design 
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The frequency response of a controller is not its only measure of suitability to the task. 
The designer must also check the time domain response to ensure adequate closed-loop 
performance. Reference [36] indicates that acceptable performance for an autopilot. is inet. 


t 


by the following conditions: 


l. No more than 5% overshoot. to a step input. 


2. 2% settling time less than 5 ship lengths at design speed. 


Figure 5-6 indicates that the autopilot design meets these criteria. Additionally, the tracking 


error at steady-state is less than 1% as dictated by the frequency response. 


Closed-Loop Step Response for Autopilot 
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Figure 5-6: Closed-Loop Autopilot Step Response 


5.2.3 Closed-Loop Simulations of the Autopilot Design 


The true test of a linear controller is how well it performs in conjunction with the nonlinear 
plant from which it was derived. Incorrect assumptions on the level of uncertainty in the 
model may lead to poor performance if the unmodeled dynamics have a strong effect on thie 
plaut. In this case. this would arise if the nonlinearities in the equations of motion have a 
strong impact on the dynamics. These nonlinearities are believed to be weak in the ship 
maneuvering problem. however, as described in Chapter 4. This proved to be true in this 


study as well. 
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To validate the controHer performance. two differeut simulations were developed. The 
first is a simulation of the controller's ability to reject. a constant current disturbance equal 
to ð kuots. For the second simulation, a maneuver was generated by computing the required 
instantaneous heading angle required for the ship to begin on a steady course and speed 
and execute a 90° left turn and steady up on the new heading. After a period of steady 


operation., the ship then executes a 90° right turn and steadies up on its new heading. 


5.2.3.1 Current Rejection Simulation 


Figure 5-7 Ulustrates the ship’s actual trajectory m relation to the dead-reckoned trajectory 
m the absence of anv disturbance. The current in this case must be simulated as an initial 
coudition on the shrp’s sway velocity, v. due to the fact that the dynamics of position. y. 


were neglected. The figure shows that the ship gets “blown” approximately 1/3 of a ship 
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Figure 5-7: Nonlinear Plant Trajectory During Current Rejection 


length off ıts initial track as a result of the current disturbance. It does. however. return 
to its conmmanded heading angle t. Thus, the controller takes the required action to return 
the ship to its intended heading. 


‘It should be noted that the standard compass rose has been rotated 90° in each of the simulations to 
follow. Thus. a heading of 000 in the figures corresponds to a heading of 090 on a standard compass rose. 
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Figure 5-8 illustrates the change 1 ship heading. The ship is initially turned to the 
right due to the side forces and moments generated by the current interaction with the ship 
hull. The controller detects the difference between the actual and commanded headings via 
noisy measurements and generates commands to the rudder to turn the ship to the left and 
return it to its intended heading. The heading angle in the figure appears to be zero. This 
is not actually the case. however. The ship actually develops a constant sway velocity that 
results ma side-shp angle. Thus. while the figure shows the ship’s actual heading equal to 
zero. the ship is actually turned at. a smal angle into the oncoming current to account for 
it. Simulations that show the actual heading angle appear in the following section where 


the dvnamies of y are included in the simulation. 
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Figure 5-8: Nonlinear Plant Heading Change During Current Rejection 
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Figure 5-9 demonstrates the commanded and actual rudder angles during the current 
rejection simulation. The commanded rudder angles may seem excessive at first glance. 
The figure shows. however. that during the maneuver, inclusion of the rudder dynamics in 
the equations of motion rejects these high-frequency oscillations. Thus, the actual rudder 
motion is smooth durimg the correction. After reaching steady-state, the rudder does exhibit. 
some high-frequency oscillatory motion. The amplitude of these oscillations. however. is less 
than 0.1° over periods of several seconds. This is not deemed to be a severe condition that. 
would tend to cause excessive wear on the mechanical components. 
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Figure 5-9: Actual and Comnianded Rudder Angle During Current Rejection 


The reason for these high-frequency oscillations mn the commanded rudder angle is illus- 
trated im Figure 5-10. Recall that the goal of loopshapmg controller design is to ensure that 
the closed-loop combination of the plant and controller attenuate high-frequency signals. 
Because the plant has very sinall gains at high frequencies, the controller is free to have 
rather large gains at these frequencies. Thus. because the input to the controller is the 
noisy. high-frequency measurement of the heading angle. the controller responds in kind. 
The controller gain at the measurement frequency is approximately 8 dB. Thus, the con- 
troller actually amplifies the measurement noise. The plant gain. however, at the same 


frequency is approximately -33 dB. Thus. the effect of the noise on the entire closed-loop 
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Sıngular Values of the Controller and Plant 
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Figure 5-10: Plant and Controller Open-Loop Gains 


system is attenuated. 


5.2.3.2 Track-Changing Maneuver 


Figure 5-11 shows the dead-reckoned and actual trajectories of the ship during a track- 
changing maneuver. The actual noulinear plant includes the effect of forward speed loss 
due to maneuvers. Thus. because the ship loses speed during the turn and experiences 
sway motion. it is not capable of completing the tight turn based solely upon kinematic 
relations. Further. the goal of the autopilot is only to maintam the desired heading rather 
than to regulate the cross-track error. A regulator of this type is described i the following 
section. The track error is approximately three ship lengths following the initial turn. This 
is reduced. however. during the second leg of the maneuver due to the sway velocity. The 
sway velocity in the, second leg is opposite that of the initial leg. Thus. the ship actually 
inoves back toward its dead-reckoned track on the second leg of the maneuver. 

Figure 5-12 shows the actual and commanded heading angle during the track-changing 
maneuver. The errors in tracking the ramp input account for the cross-track errors in the 


trajectory. The tracking error in steady-state. however. is again less than 1%. 
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Figure 5-11: Nonlinear Plant Track-Changing Maneuver with Autopilot 
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Figure 5-12: Noulinear Plant Heading Angle during Track-Changing Maneuver with Au- 
topilot 
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Figure 9-13 shows tlie commanded and actual rudder angles during tle track-changing 
maneuver. Again. the rudder motion is quite smooth during the major portions of the 
maneuver and is not affected by the high-frequency controller commands. Further, notice 
that the rudder sweep during tle maneuver is approximately 20° over a period of about 
three minutes. Thus. the controller produces a slow, smooth rudder motion for the duration 
of the maneuver. The steady-state oscillations are again small and not deemed to cause 
excessive mechanical wear. 
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Figure 5-13: Actual and Commanded Rudder Angle During Track-Changing Maneuver 
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5.3 Cross-Track Error Controller Design 


The goal of the cross-track error controller is to regulate the ship’s lateral deviation from 
its intended track. The previous section demonstrated that an autopilot is not sufficient for 


this purpose. The controller must meet all criteria previously stated. 


5.3.1 Dynamic Model 


To begin tlhe design of the cross-track error controller, tlle coordinate system was shifted. 
In the new coordinate system, the x-coordinate represents tlie distance traveled along the 
reference trajectory aud the y-coordinate represents tlie lateral deviation from the intended 
track. The heading angle now represeuts the deviation from the desired heading angle along 


the intended track. The new coordinate system is illustrated in Figure 5-14. The kinematic 
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Figure 5-14: New Coordinate System for Track Controller 


relations governing ship motions in this new coordinate system is given by the following 


equations: 


Tt = u cos yY — vsin y (00) 


y = usinw + vcos w (5-30) 
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Linearizing aud non-dimensionalizing these equations leads to equation 5.37. 


y= (5.37) 


where the dynamics of sway have again been removed as in 5.2.1. The dynamics of x have 
also been removed since it is not a control variable. Thus, differentiating equation 5.37 
aud inserting the result in equation 5.31 leads to the following expression for the control 


variable. y: 
YT 1038) 


Thus. using SISO loopshaping control design theory, the track controller may be designed by 
sunply adding an integrator to the plant transfer function given by equation 5.29. Thus. the 
new plant transfer fiction that maps the control input, ðe. to the lateral track deviation. 


y. becomes: 


P 
om (5.39) 


Py = 


5.3.2 Loopshaping Controller Design 


Plants with double integrators can be problematic for the control systems engineer due to 
the phase shift introduced by the two poles at the origin. The two poles introduce a 1809 
phase shift in the plant at DC. Thus. because the dynaniics of yaw rate. r, behave as a 
second-order system. the phase decreases as frequency increases because the plant has słow 
poles. 

The slope of the loopshape curve as it passes through the crossover frequency (where 


the magnitude is one) must satisfy the condition given by equation 5.40. 


dL 1 
E er. (5.40) 
dw a 


The phase associated with tlis condition is -180°. System stability requires that the phase 
be greater than this when the magnitude is one. Thus, the closed-loop system will be 
unstable if this condition is not met. Thus, the simple loopshaping design described in 5.2.2 


may uot be employed. In this instance. the controller must have some zeros to add the 
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required phase around the crossover frequency. Thus. the controller design must proceed 
by direct design of the controller. C *, computation of S and T, and verification of the 
desired frequency respouse of the loop gain. L = P,C. Thus, the general controller transfer 


function to start the design may be defined as follows: 


_ k(s + 21)(s + 29) er (s + Zn) 


Sal 
(< Fanals + po)... (s + Pn) on 


where nm may be as many poles and zeros as required to achieve a desirable loopshape, 
subject to the requirement that the transfer function remain proper (degree of numerator < 
degree of denominator) [11]. The final controller design may be described by the following 
transfer function: 


20(s + 0.5)%(s + 6)(s + 1000) 


5 ETS + 15)(s + 30) 


(5.42) 


Figure 3-15 illustrates the loopshape for the cross-track error controller design. The 
design proceeded iteratively by adding zeros and poles im the controller to achieve the 
desired frequency response properties. Figure 5-16 additionally shows the gain aud phase 
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Figure 5-15: Loopshape for Cross-Track Error Controller Design 


“Zeros and poles may be easily added and modified using the zpk function available in MATLAB. 
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wargins assoclated with the design. Again. the margius fall well witliin the range of accepted 
practice. Notice im the phase plot of the Bode diagram how the controller adds phase to 
the svsten at frequencies below the crossover frequency, thus stabilzing the plant. 
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Figure 3-16: Gain and Phase Margins for Cross-Track Error Controller Design 





Figure 0-17 illustrates the system closed-loop step response. The overshoot and settling 
tine in This case may seem excessive when compared to conventional design practice. Com- 
pare. however. Figures 5-17 and 5-6. and consider equation 5.37. It would be extremely 
difficult. if not impossible. to command the dynamics of y to be faster than tlre dynamics 
of y because y depends solely upon y. Thus, its dynamics must be inherently slower. This 
plivsical fact is confirmed by the step response plot. Closed-loop simulation will subse- 
quently show. however. that this controller performs extremely well on the nonlinear plant 


with smooth maneuveriig trajectory commands. 
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Figure 5-17: Cross-Track Error Controller Closed-Loop Step Response 
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5.3.3 Closed-Loop Simulation of the Cross-Track Error Controller 


By shifting the coordinate system to the one shown in Figure 5-14, the two simulations 
performed in 5.2.3 may be combined into one simulation. The simulation speed is again 22 


knots as in the previous simulations. The simulation steps are as follows: 


1. The ship begins on a steady course and speed with no current. 


o 


After a short time, the ship experiences a current gust of 5 knots directly on its port 


bean. 


3. Shortly thereafter. the slip makes a left turn, moves two ship lengths to port and 


steadies up. 
4. After a period of steady travel. the ship executes a right turn and returns to its original 
track. 
5.3.3.1 Maneuvering Simulation in Calm Seas with Current 


Figure 5-18 shows the desired and actual trajectories as well as the heading angle during the 
maneuver. At the beginning of the maneuver, the controller takes strong action to reject 
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Figure 5-18: Trajectory and Heading during Nonlinear Plant Simulation 


the current disturbauce. It. overshoots the desired track slightly upon its return and begins 
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to compensate. The left turn then executes and the ship follows the desired track with very 
small error until it steadies up. Notice that during the steady parts of the maneuver. the 
slup manitams a constant heading offset to counteract the current. 

Figure 5-19 shows the commanded and actual rudder angle during the maneuver. The 
controller is again seen to respond significantly to the sensor noise for the same reason listed 


i 5.2.3.1. The most noteworthy item in this figure. however, is the way the controller senses 
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Figure 5-19: Actual and Commanded Rudder Angle During Nonlinear Plant Simulation 


the current aud uses it to execute the right turn to return to its original track. Note how 
the ship's heading decreases slightly below its offset. but does not pass through zero. Thus. 
the controller uses the current to allow the ship to be “blown” back onto its desired track. 

Some may note that the ship requires a great deal of time (approximately 8 minutes) to 
complete the turning maneuvers. This may or may not be desirable. Smooth trajectories are 
required. however, to prevent the controller from over-compensating during the overshoot 
phase. This application could be extremely valuable. however, in maneuvers to avoid other 
vessels. For example. modern radar is capable of identifying contacts well before (longer 
than 8 minutes) thei closest point of approach (CPA) time. Maneuvers to avoid other 
vessels are currently calculated by watchstanders using vector diagrams. Thus. a digital 


computer could calculate the CPA, compute a suitable track to avoid the contact by some 
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threshold distance. and then compute a snivoth trajectory to maneuver to the new track. 


observe CPA. and maneuver back to the intended track. 


5.3.3.2 Maneuver in Sea State 4 with 5-Knot Current 


The simulation would not be complete without verifying the controller’s performance in a 
sea state. The sea state is modeled as a random disturbance and filtered to preserve only the 
low-frequency components. Table 6 of [38] indicates that a state four sea has a significant 
wave height of 1.25-2.5 meters (4 - 8 feet). Additionally, Table 3 of [38] indicates tliat waves 
of this height represent sea conditions with the highest probability of occurrence worldwide. 
This. sea state four was chosen for simulation Figure 5-20 shows the wave spectrum and 


Figure 5-21 shows the time history of the wave heights. 
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Figure 5-20: Power Spectrum of Wave Model 
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Figure 5-21: Time History of Wave Heights 
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The same maneuver described in 5.3.3 ıs performed with the current aud sea state 
introduced as disturbances. Figure 5-22 slıows the actual and commanded trajectories as 
well as the heading angle during the nianeuver. Notice that. smmiar to the maneuver 11 
calin seas. the controller maintains a constant heading offset to counteract the effects of the 
current. The low-frequency wave disturbances. however, cause the ship to deviate from its 
connuauded track. Thus, the ship oscillates a small amount about its commanded track. 


The amplitude of the oscillations, however. is less than 20 yards (about 15% of the ship 
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Figure 5-22: Trajectory and Heading During Track Change 
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Figure 5-23 shows the commanded and actual rudder angles during tlie nıaneuver. Care- 
ful examination of Figure 5-23 mdicates a low-frequency signal among the ligh-frequenucy 
measurements driving the controller. It is this low-frequency command signal that allows 


the ship to remain on track despite the disturbances. 
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Figure 5-23: Commanded and Actual Rudder Angle during Track Change 
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This low-frequency component may be seen more clearly in Figure 5-24. This figure 
shows the time history during the steady portion of the maneuver between the two turns. 
Notice how the rudder corrective actions exhibit a period simular to that of the ship devia- 
tious about its intended track. In particular, the figure shows the rudder’s large corrective 
action applied at the point about 1400 seconds into the simulation. This is the point where 
the ship exhibits its largest track deviation due to overshoot at the end of the turn. The 


controller seuses this. however. and takes the appropriate action. 
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Figure 5-24: Low-Frequency Component of Actual Rudder Angles 
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Chapter 6 


Conclusion 


61 Summary 


This work gives rise to a systematic method for adding automatic maneuvering coutrol 
systein desigu to the macroscopic slip design process. The procedure may be summmarized 


as follows: 


Build a scale model of tle ship. 


—" 


to 


Outht the model witli seusors. 
3. Conduct a series of controlled maneuvers with the model to gather data. 


4. Design a SAEKF and use the data collected from the model to identify the ship 


hydrodynamic coefficients. 
5. Use the identified coefficients to develop controllers. 
6. Test the controller designs on the scale model. 


This work accomplishes a piece of this process by identifying the linear damping derivatives 
and developing control laws based upon them. Specifically, this work accomplished the 
following objectives: 

e Selection of an assumed form of the equations of ship motion in the horizontal plane. 

e Design Of a SAEKF to identify the linear damping hydrodynamic coefficients within these 
equations. 


e Successful identification of these coefficients which are the most critical to overall ship 
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maneuvering stability in the horizontal plane. 
e Design of two simple controllers that may be tested on scale ship models to determine 
thei performance. 

The current method of identifying these terms involves large towimg tanks and rotating 
arm basins. It also requires extensive labor to operate the equipment. These requirements 
are in addition to the need for a scale model of the ship. Thus, development of the ideas 
mtroduced in [3] and [43] and continued here may not only lead to an efficient, effective 
process for adding maneuvering automation to the ship design process, it may also reduce 


the infrastructure required to support it. 


6.2 Conclusions 


The following list summarizes the items deemed most noteworthy by the autlıor over tlıe 


course of this stndy. 


l. The EKF method shows definite feasibility as a method to determine linear slip 
hydrodynamic coefficients. These parameters may then form the basis for simple 


coutroller designs that can be tested on scale models. 


2. The parameter estimates provided by the SAEKF need not be very accnrate to pro- 
vide values suitable for controller design. The parameters determined in this study 
reproduced simulation trajectories with approximately 15% error in steady turning 
diameter. Errors as large as 40%, however. may very likely yield acceptable controller 


designs for slow maneuvers. 


3. The dynamics of sway need not be considered to develop effective linear controllers 
for the types of manenvers considered in this study. Thus, controller designs for these 


types of maneuvers may be based on only two parameters. 


4. The EKF performs extremely well as a state estimator (due to the plıysical state 
measurements) despite significant inaccuracies in the filter dynamic model. Specifi- 
cally. the SAEKF model used in this study was completely linear, but still produced 
very accurate estimates of the physical states (u, v, r, and n) of the nonlinear plant 
model used to generate the measurement data. The dynamics of u and n, in fact., 


were completely ignored in propagating the filter state estimates. 
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9. The perforinance of the SAEKF as a parameter estimator, however, will degrade if the 


filter model is not sufficiently accurate in describing the actual plant. The modeling 


errors will cause significant bias in the parameter estimates. 


The instantaneous rudder angle of attack decreases rapidly with respect to the nominal 
augle of attack as sway speed and yaw rate increase. Thus. the linear lift coefficient 
assunption proved to be valid in the nianeuvers conducted in this study despite the 


use of rudder angles traditionally believed to lie outside the linear regime. 


6.3 Recommendations for Further Study 


As in any other design study. this work leaves several areas that require further thought and 


development. Among these. the most important in the author’s mind may be summarized 


as follows: 


E 


IN 


Gi 


This work considers performing the identification in calm water. Environmental dis- 
turbances such as wind and waves could have a significant effect on the identification 


process. The magnitude of these effects should be investigated. 


Rolling motions have an effect on maneuvering in the horizontal plane. Thus. the 
inclusion of rolling motions in the nonlinear model could have an effect on the linear 


controller performance. The maguitude of these effects should be investigated. 


Despite the identification of four terms in this study. these four terms were reduced to 
two to develop the controllers. Identification of only these two terms using an SAEKF 


may be possible. This could make the implementation of the SAEKF mucli simpler. 


The EKF provided very accurate estimates of the surge speed and propeller rotational 
speed despite the fact that its model completely ignored their dynamics. These quai- 
tities are easily measured by the speed log and a shaft tachometer. Thus, perhaps the 
EKF could be used to provide state estimates for a MIMO controller to control the 


slıip's x-coordinate as well as its y-coordinate. 


The proper noise covariance terms are an important design parameter in any Kalman 
filter implementation. This study determined the correct noise covariances through an 


iterative procedure of identification and simulation. Inclusion of the noise covariances 
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as filter parameters to be estimated may reduce the iterative nature of this process 


aud lead to a more efficient design process. 


linplement the method outlined herein on the DDG-51 scale model developed at the 


MIT Towing Tank. 


hhivestigate the range over whicli the assumption of linearized lift on rudders is valid 


by accounting for the ship dynamics in computing the instantaneous angle of attack. 


00 
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Appendix A 


State Augmentation of the 
Extended Kalman Filter 


The derivations presented in section 2.2 hold for estimation of the state vector describ- 
ing a nonlinear svstem. Auginentation of the state vector for use in estimating unknown 
paraieters is the subject of the sequel. 

Recall equations 2.10 and 2.19 in 2.2.2 and 2.2.3, respectively. The variable x represents 
the estimate of the states describing the physical system. Therefore, the estimated state 
vector can be augmented with the unknown parameters im the system to allow then to be 
estimated along with the system states. In the case of DDG-51, the unkown parameters are 
coustant hydrodynamic maneuvering coefficients. Maneuvering coefficients may change with 
changes iu subinerged geometry. speed. etc. Nondimensionalization of the equations removes 
speed considerations and submerged geometry certainly changes slowly with time. DDG- 
ol’s compensated fuel systen1. however. ensures that the ship maintains nearly constant 
draft over its entire range of operation. Therefore, the state vector is augmented with 
the nuknown system parameters whose derivatives are zero. For example, to estimate ax 


uuknowu parameter. 6. in a system fully described by states 11 and xo, the state vector amd 
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Appendix A 


State Augmentation of the 
Extended Kalman Filter 


The derivations presented in section 2.2 hold for estimation of the state vector describ- 
ing a nonlinear system. Augmentation of the state vector for use in estimating unknown 
parameters is the subject of the sequel. 

Recall equations 2.10 and 2.19 in 2.2.2 and 2.2.3, respectively. The variable £ represents 
the estimate of the states describing the physical system. Therefore, the estimated state 
vector can be augmented with the unknown parameters in the system to allow them to be 
estimated along with the system states. In the case of DDG-51, the unkown parameters are 
constant hydrodynamic maneuvering coefficients. Maneuvering coefficients may change with 
changes in submerged geometry. speed. etc. Nondimensionalization of the equations removes 
speed considerations and submerged geometry certainly changes slowly with time. DDG- 
ols compensated fuel system, however. ensures that the ship maintains nearly constant 
draft over its eutire rauge of operation. Therefore, the state vector is augmented with 
the unknown system parameters whose derivatives are zero. For example, to estimate an 


unknown parameter. b. in a system fully described by states zı and x», the state vector and 
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associated F matrix would be augmented as follows: 


1 ENO t) 





Load) gel) 

b 0 

| oh ay a2 
0 0 


Using this methodology. an initial estimate is made for the unknown parameters in the 
augmented staic vector aud the associated error iu the estimate. Thi parameters are then 
estimated iu exactly the same manner as the state vector. Once the state vector has been 
angmented with the unknown parameters. and the F matrix modified appropriately. tle 


EKF becomes the SAEKF. The following section presents an example application of the 


EKF. 
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A.l Application of the State Augmented Extended Kalman 


Filter to a Nonlinear Tracking Problem 


The Extended Kalman Filter has proven to be extremely effective for tracking applica- 
tions (such as tracking an object. with a radar). The radar provides position information 
for the object. but no velocity information. Effective tracking. however. requires velocity 
formation to fully describe the state of the system. 

Assume for suuplicity that an object falls through space directly toward a tracking radar 
as illustrated in Figure A-1 (i.e. a one-dimensional tracking problem). Assume furtlier that 
the system has an unknown parameter (tle ballistic coefficient. b) which affects the drag 
force on the object as it falls through space. The goal is to estimate the velocity of the 
object. v. and the unknown parameter, b, with only noisy measurements of its position. z£. 
available as actual data. The nonlinearity arises through the dependence of drag on the 


square of the object's velocity and the ballistic coefficient (d(b) x v7). 
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Figure A-1: A One-Dimensional Nonlinear Tracking Problem 


‘This example is taken from example 6.1-2 in [19]. 
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Figure A-2 illustrates the results. The upper plot shows the filter’s ability to very 
accurately estimate the object's velocity as it falls through the atmosphere. The lower plot 
ustrates the filters ability to converge ou the true value of the ballistic coefficient using the 
measurements aud computed gains. The filter does not change much iitially because the 
object experiences very little drag high in the atmosphere. As the object enters the dense 
atmosphere the filter converges to the true value and remains there. Thus. one concludes 
that the filter has difficulty identifying the parameter value if it does not significantly affect 


the dvnannes. 
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Figure A-2: EKF Application 
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